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Abstract . 


After  a  critical  survey  of  existing  theories  of  erosive  burning 
and  of  published  experimental  data,  the  combustion  mechanism  of  ammonium 
perchlorate  propellants  is  discussed  and  a  model  of  erosive  burning  of 
these  formulations  is  proposed. 

Erosion  is  assumed  to  smooth  the  surface  whose  burning  rate  is 
then  controlled  by  oxidizer  regression,  no  energy  being  supplied  by  the 
main  diffusion  flame.  Furthermore  the  adiabatic  decomposition  temperature 
of  the  perchlorate  is  assumed  to  be  obtained  at  the  end  of  the  laminar 
sublayer  whose  thickness  can  be  related  to  mean  flow  velocity  through  fluid 
dynamics.  Correlation  of  test  data  permits  to  determine  the  transition 
point.  Grain  design  and  motor  scale  are  taken  into  account  through  local 
hydraulics  radius. 

Low  pressures,  large  radii  and  high  burning  rates  are  shown  to 
reduce  erosive  burning. 

Finally  an  attempt  is  made  to  explain  irregular  burning  by- surface 
nature  and  roughness,  through  erosivity. 
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1 .  INTRODUCTION. 


Experimental  data  show  that  the  "burning  rate  of.  a  solid  propellant 
is  increased  above  its  normal  pressure-dependent  value  by  the  presence 
of  a  flow  parallel  to  the  burning  surface. 

This  phenomenon,  which  plays  an  important  role  in  internal  bal¬ 
listics  of  high  performance  rockets,  is  known  as  EROSIVE  BURNING. 

Its  existence  has  been  recognized  early  in  the  developmentt  of 
solid  propellent  motors.  Despite  several  important  contributions, 
however,  it  has  not  deserved  much  attention  until  recently,  since  it 
has  been  possible  to  overcome  rather  easily  by  cut  and  try.  its  major 
detrimental  effects  on  relatively  small-sized  motors. 

The  cost  of  firing  the  very  large  solid  propellant  rockets 
presently  in  the  development  phase  precludes  this  empirical  approach. 

Moreover  the  method  permits  to  overcome  only  the  most  obvious 
disadvantages  resulting  from  erosive  burning,  such  as  the  appearance 
of  a  pressure  peak  after  ignition,,  or  burnthrough  due  to  uneven  combustion 
along  the  grain  length  and  premature  contact  of  the  hot  gas  with  the 
chamber  wall  at  the  nozzle  entrance 'where  the'  flow  velocity  is  a 
maximum. 


The  most  severe  difficulty  in  studying  these  phenomena  is  the 
absence  of  a  sufficient  volume -of  .reliable  and  systematic  experimental 
information  wrhich  results  in  representing  erosive  burning  by  means  of 
oversimplified  relationsnips  which  barely  fit  the  data  for  one  propellant 
but  do  not  even  describe  the  behaviour  of  another  similar  formulation. 
These  empirical  relationships,  which  assume  linearity  between  burning 
rate  increase  and  flow  velocity,  Mach  number,  mass  velocity  or  reduced 
mass  velocity,  may  rather  satisfactorily  represent  the  phenomenon  under 
severe  erosive  conditions.  In  moBt  cases,  -however,  they  appear  to  be 
insufficient  to  predict  the  burning  rate  u^dsr  conditions  of  lean 
erosion,  even  if  the  existence  of  a  controversial  "threshold  velocity" 
is  taken  into  account. 

In  fact  the  analysis  of  available  data  suggests  that  the  burning 
rate  dependence  on  flew  velocity,  Mach  number,  mass  velocity  or  reduced 
mass  velocity  can  exhibit  many  patterns  ,  depending  on  propellant 
formulation  and  somewhat  on  combustive  pressure.  Figure  1  represents 
four  possibilities  which  can  be  encountered. 


FLOW  VELOCITY  .  MACH  NUMBER  , 
MASS  VELOCITY  or  REDUCED  MASS 

VELOCITY 
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Closer  analysis  shows  that  the  burning  rate  without  erosion  '  appears 

more  or ’ less  as  an  anomaly.  Under  strong  erosion  it  is  doubtful  that  the 
parameters  a  and  n  of  the  classical  law: 


still  play  a  predominant  role.  " 

Any  empirical  relationship  giving  the  total  burning  rate  r,  inclu¬ 
ding  erosion  effect,  in  terms  of  r  is,  therefore,  suspicious,  even  if 
its  use  appears  to.  provide  valuable  results.  This  point  of  view  has 
been  developed  in  a  previous  report  dealing  with  the  influence  of  erosive 
burning  on  internal  ballistics  f  i  J 


Prediction  of  the  influence  --of  such  factors  as  grain  geometry  and 
scale  or  minor  changes  in  formulation  on  erosive  burning  could  only  be 
achieved  through  a  fair  physical  understanding  of  the  phenomenon,  taking 
into  account  the  combustion  mechanism  of  the  propellant. 

Moreover  an  improved  understanding  of  the  burning  mechanism  and  of 
the  erosivity  may  provide  a  tentative  explanation  fop  unstable  combustion 
in  solid  propellant  rockets. 

Despite  years  of  intensive  research  by  many  groups,  these  problems 
are  still  barely  understood  :  the  speculative  nature  of  the  present  work 
is  thus  plainly  realized  by  the  author.  '  Indeed,  over  and  above  several 
drastic  assumptions,  the  theory  developed  hereunder  is  entirely  based  on 
a  simplified  combustion  mechanism  on  which  a  rather  rough  fluid  mechanics 
model  has  been  superposed  for  taking  -into  account  erosive  burning.  The 
combination  of  these  two  theories,  however,  provides  some  interesting 
possibilities  for  cross-checking  the  assumption  with  the  experimental 
results;  by  all  means  this  type  of  frontal  approach  is  the  only  possible 
one  for  progressing  towards  the  understanding  of  the  phenomena. 


XXX 


In  this  report,  special  emphasis  has  been  given  to  the  case  of  composite 
propellants  using  ammonium  perchlorate  as  oxidizer,  in  view  of  their  wide 
range  of  application. 


2.  SURVEY  OF  EXISTING  THEORIES  OF  EROSIVE  BURNING. 


This  survey  does  not  intend  to  describe  in  details  the  mathematics 
of  the  various  theories  which  have  been  proposed  for  explaining  erosive 
burning,  but  rather  to  discuss  the  hypothesises  on  which  these  models 
are  based. 

All  these  theories  neglect  the  mechanical  effects  of  flow  on  the 
propellant  and  assume  that  the  increase  in  burning  rate  results  from 
accelerated  chemical  reactions  under  the  influence  of  some  mechanism 
of  fluid  dynamics.  These  theories  can  be  divided  in  two  broad  categories 
depending  on  whether  the  flame  structure  is  taken  into  account  or  not. 


A.  The  Lenoir  -  Robi'llard  Theory. 


The  best  known  theory  of  erosive  burning  of  solid  propellants  has 
been  proposed  by  Lenoir  and  Robillard  [  2  ]  .  It  has  also  been  surveyed 
by  Green,  Penner  and  Schultz  who  described  it  as  a  major  advance  towards 
a  description  of  the  erosive  burning  mechanism.  [31  . 

Lenoir  and  Robillard  have  applied  an  empirical  relationship  of  the 
type  employed  by  Rannie  and  others  in  studies  of  heat  transfer  to  trans¬ 
piration-cooled  surfaces,  under  the  assumption  that  the  presence  of  a 
flow  parallel  to  the  surface  accelerates  the  convective  heat  trani  er 
between  the  main  gas  core  and  the  burning  surface . 

It  is  believed  that  this  model  suffers  from  a’ fundamental  drawback: 
the  basic  ■•'elationship  used  by  Lenoir  and  Robillard  describes  the  heat 
transfer  to  transpiration  cooled  surfaces  only  when  the  temperature 
gradients,  in  the  gas.  are  due . t o  convection  only,  while  in  the  presence 
of  a  burning  surface  the  temperature  gradients  are  determined  primarily  by 
the  chemical  reactions  in  the  combustion  zone,  the  flow  pattern  having 
only  a  second  order  influence. 

The  good  numerical  results  are  not  an  absolute  proof  of  the  validity 
of  their  model  since  two  constants  are  computed  from  experimental  data  in 
order  to  obtain  the  agreement. 

This  theoretical  objection  is  supported  by  the  observations  of  several 
authors  who  reported  that  erosive  burning  is  not  influenced  by  the  temper¬ 
ature  of  the  main  gas  flow  |  4,  6,  7,  14]  ,  and  even,  ir.  certain  cases 
under  sufficiently  strong  erosive  conditions,  by  the  combustion  temperature 
of  the  propellant  itself  (14]- 

The  theory  proposed  by  Lenoir  and  Robillard  has  been  used  successfully 
for  explaining  such  behaviours  as  the  inverse  dependence  of  the  erosive 
effect  upon  propellant  burning  rate  or  the  controversial  existence  of  an 
apparent  velocity  threshold  near  the  fore  end  of  the  grain. 


We  believe,  however,  that  this  success  is  due  to  the  fact  that  only 
a  few  cases  of  practical  importance  have  been  studied  without  attempting  a 
systematic  investigation  of  the  phenomenon  in  the  full  range  of  pressure, 
formulation,  granulation  or  scale. 


B.  The  Corner  Theory  and  its  Development  a . 


As  pointed  out  above,  the  burning  rate  is  primarily  determined  by 
chemical  reactions  and  erosion  introduces  a  second  order  effect  only. 

A  theory  of  erosive  burning  must,  therefore,  take  into  account,  the  flame 
structure  and  describe  how  it  is  affected  by  fluid  dynamic:').  Such  a 
frontal  approach  is,  basically,  difficult  since  it  involves  the  coupling 
between  two  phenomena  which  ,  even  separately,  are  barely  understood 
i.e.  the  combustion  mechanism  and  the  flow  field. 

o  Corner. 


The  first  approach  along  these  lines  is  due  to  Corner  [41. 

After  having  remarked  that  near  the  aft  end  of  the  charge ,underr  ^ 

severe  erosion,  typical  values  of  the  Reynolds  number  range  from  ICr  to  10  , 
this  author  used  the  Prandtl-Karman  momentum  theory,  for  representing. the 
boundary  layer  and  for  computing,  within  the  combustion  zone,  the  increase- 
In  thermal  conductivity -which  is  introduced  by  turbulence,  above  the 
laminar  sublayer.  5 

This  appr.oaoh  gives  the  effective  thermal  conductivity  which  ■  4:es 
into  account  eddy  diffusiviiy,  as  a  function  of  the  distance  from  the  ■ 
burning  surface  as  well  as  the  mean  flow  velocity  in  the  section  under 
consideration. 

Therefore  it  should  be  possible,  in  principle,  to  integrate  the 
equations  describing  the  combustion,  with  a  variable’  conductivity,  at  the 
price  of  a  considerable  increase  in  complexity. 

Alternatively  a  suitable  mean  value  of  the  distance  from  the  surface 
night  be  selected  that  would  give  a  mean  conductivity  applicable  in  the 
most'  important  region  of  the  flame.  Corner  selected  the  distance  at  which 
the  reaction  is-half  completed  without  erosion,  an  assumption  which  does 
not  take  into  account  the  fact  that  this  characteristic  distance  almost 
certainly  decreases  a c  the  mean  flow  velocity  increases. 

Altough  no  serious  attempt  was  made  to  draw  systematically  quantitative 
results,  Corner's  theory  provides  a  good  explanation  for  the  following 
experimental  facts: 

-  slow  burning  propellants  are  more  sensitive  to  erosion  than  fast  burning 
ones  which  are  characterized  by  thinner  flames, 

-  erosive  burning  is  insensitive  to  the  temperature  of  the  main  flow  since 
the  surface  is  bathed  in  its  own  products  and  these  are  set  into  trans¬ 
verse  motion  by  the  main  str 


earn  , 
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-  the  threshold  velocity  ‘below  which  erosion  does  not  increase  the  burn¬ 
ing  rate  is  due  to  the  penetration  of  turbulence  in  the  combustion  zone. 

Two  major  objections,  however,  can  be  drawn  against  Comer's  approach. 

First,  from  the  point  of  view  of  fluid  mechanics,  even  its  limited 
success  appears  remarkable  since  it  is  based  upon  empirical  expressions 
for  the  velocity  distribution  in  a  one-dimensional  flow  adjacent  to  an 
inert  boundary  and  the  effects  of  mass  addition  normal  to  the  stream  are 
not  considered.  This  type  of  one-dimensional  approximation  would  be 
expected  to  yield  a  reasonably  accurate  representation  of  the  actual  flow 
field  only  near  the  aft  end  of  the  grain  where  the  ratio  of  normal  injection 
velocity  to  main-stream  velocity  is  lowT,  of  the  order  of  0.01  r.  3  i  . 

» 

Furthermore  the  one-dimensional  velocity  distribution  is  valid  only 
under  the  assumption  of  a  fully  developed  turbulence  for  which  the  char¬ 
acteristic  dimension  in  the  Reynolds  number  is  the  pipe  diameter.  In 
internal  ballistics  the  length  to  diameter  ratio  of  the  perforation  is 
seldom  sufficiently  large  for  neglecting  the  two-dimensional  nature  of 
the  flow  field  which  is  influenced  by  the  distance . from  the  fore  end  of 
the  grain . 

Secondly  Corner  assumed  that  the  rate  controling  step  of  the  com¬ 
bustion  is  the  main  gas  phase  reaction.  This  may  he  a  realistic  assumtion 
at  the  very  high  pressures  encountered  in  guns,  ahove  1.000  kg/cm2. - 
In  the  normal  pressure  range  used  in  rocket  motors,  below  100  kg/cm  ,  it 
is  generally  agreed  that  the  rate  controlling  step  is  the  solid-phase 
decomposition  governed  by  the  combustion  surface  temperature  T  through  an 
Arrhenius  type  relationship. 

Furthermore  the  structure  of  the  flame  is  rather  complex  and  in  many 
cases  several  zones  can  he  distinguished.  ’ 

P  -  Boisson  and  Tavernier. 


In  an  attempt  to  correlate  the  experimental  data  obtained  by  Berger, 
Prache  and  Tavernier  t  5  i  >  Boisson  and  Tavernier  extended  the  Corner 
theory  for  taking  into  account  the  mass-diffusion  due  to  turbulence,  in 
addition  to  the  increase  in  thermal  conductivity  t  6 ,  7  1  •  These  authors 
concluded  that  in  the  low  pressure  range,  the  turbulent  mass-diffusion  has 
a  predominant  influence  on  erosive  burning.  No  attempt  has  been  made,  how¬ 
ever,  for  systematically  exploring  the  quantitative  influence  of  the  various 
parameters . 


Although  Boisson  and  Tavernier  contributed  to  the  physical  understanding 
of  erosive  turning  by  emphasizing  the  importance  of  turbulent  mass-diffusion 
in  the  reaction  zone,  their  theory  suffers  from  the  same  drawbacks  as  the 
Corner  theory:  the  fluid  mechanics  model  is  oversimplified  and,  moreover, 
the  assumption  of  a  main  gas  reaction  controling  the  burning  rate  does  not 
appear  realistic  in  the  low  pressure  range. 
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y  -  Geckler. 


first  serious  attempt  to  take  into  account  the  flame  structure  for 
computing  the  erosive  burning  constant  is  due  to  Geckler.  [8  j  . 

In  his  famous  survey  about  combustion  mechanisms  of  solid  propellants 
this  author  adapted  the  Corner  theory  for  taking  into  account  the  flame 
structure  cf  double  base  propellants  which  has  been  developed  by  Parr,  Craw¬ 
ford,  Rice  and  C-insll.  According  to  this  model,  the  rate  controling’  step 
of  the  corr.Dust.ion  is  the  solid— phase  decomposition  which  is  governed  by  an 
Arrhenius-type  equation,  and  the  gas-phase  reaction  can  be  divided  in  three 
parts.  Close  to  the  surface  exothermic  reactions  take  place  in  the  fizz 
sche.  Then,  in  the  preparation  zone,  activated  products  are  formed  without 
heat  production.  Finally,  when  a  sufficient  concentration  of  activated 
products  is  achieved,  the  final  reaction  occurs  in  the  combustion  zone  at 
the  end  of  which  the  isobaric  combustion  temperature  is  obtained  (figure  2). 
Furthermore  these  authors  have  shown  that  at  sufficiently  low  pressure, 
below  100  atm,  the  thickness  of  the  preparation  zone  is  mucli  greater  than 
the  thickness  of  the  ether  zones,  and  the  heat  transfered  from  the  flame 
zone  into  the  fizz  zone  can  be  neglected.  Geckler  thus  selected  the  thick¬ 
ness  y ^  of  the  fizz  zone  as  the  suitable  mean  value  of  the  distance  from 
the  surface  representing  the  most  important  region  of  the  flame. 

Using  for  y^  the  values  obtained  by  Rice  and  Ginell,  Geckler  computed 
for  that  particular  distance,  through  the  Prandtl— ICarman  momentum  theory, 
the  effective  thermal  conductivity  taking  into  account  the  turbulent  term 
and  applied  the  result  to  Corner's  equation: 

*  I 

r _  _  /  _a_  \  o .  5 

r  o  1  * 

with  r  the  actual  rate  under  erosive  conditions, 

a  the  corresponding  burning  rate  without  erosion 
A  the  thermal  conductivity  without  turbulence , 
and  a '  the  effective  thermal  conductivity  at  distance  y  ^  taking  into 
account  turbulence. 

In  applying  this  method  to  propellant  HES  4016  which  has  been  studied 
in  great  details  by  Rice  and  Ginell,  Geckler  obtained  an  extraordinary  good 
agreement  with  the  experimental  data.  He  remarked,  however,  himself  that 
although  remarkable,  this  agreement  is  probably  fortuitous.  Indeed,  from 
the  fluid  mechanics  point  of  view,  this  treatment  suffers  from  the  same 
drawbacks  as  Corner's  theory  discussed  previously.  Furthermore ,  the  preceding 
relationship  has  been  established  by  Corner  in  the  assumption  of  a.  gas-phase 
reaction  controlling  the  burning  rate.  This  is  not  consistent  with 'the 
assumption  of  the  solid-phase  decomposition  control  of  the  burning  rate  which 
has  been  emphasized  through  Geckler 's  survey.  Nevertheless  this  author 
remarked  that  this  equation  provides  a  reasonable  approximation.  For  check¬ 
ing  its  validity,  we  made  a  rough  calculation  which  seems  to  indicate  that 
the  exponent  0.7  may  be  mere  realistic  than  0.5. 


TEMPERATURE 


FIG. 2 


Finally,  Geckler  remarked,  that  this  theory  suffers  from  the  draw¬ 
back  of  not  ts.king  into  account  the  decrease  in  the  fizz  zone  thickness 
which  results  from  an  increase  in  burning  rate  due  to  erosion. 

S  -  Vandenkerekhove . 


In  attempting  to  improve  the  theory  further,  Vandenkerekhove 
proposed  a  complex  model  of  erosive  burning  based  on  rather  controver- 
sial  phenomenological  assumptions,  which  give  fair  quantitative  results 
[  9'l  •  After  having  established  the  general  heat  transfer  equation  across 
the  fizz  zone  of  a  double  base  propellant,  for  th*  purpose  of  relating 
the  surface  temperature  T  to  the  zone  thickness  y.  and  to  the  temperature 
T^  reached  at  the  end  of  the  zone,  he  remarked  that  this  relationship  is 
valid  in  the  laminar  boundary  layer,  even  in  the  presence  of  a' strong 
velocity  gradient  parallel  to  the  surface,  but  that  in  the  turbulent 
region  great  caution  is  necessary. 

Nevertheless  one  can  expect  that  the  flame  structure  still  holds 
in  the  transition  sublayer  where  turbulence  is  not  yet  fully  developed 
and  whero  tho  dimensions  of  tho  oddioo  nro  much  omallor  than  the  zone 
thickness. 

It  seems  therefore  reasonable  to  assume  that,  at  least  in  the  low 
velocity  range,  temperature  is  not  signif icantly  influenced  by  turbulen¬ 
ce  in  the  preparation  zone  across  which  the  heat  transfer  remains  negli¬ 
gible  . 


Furthermore  the  influence  of  turbulence  is  ..ot  only  to  increase  the 
thermal  conductivity  by  transporting  hot  gases  towards  the  surface. 

Hitherto  unreacted  products  are  also  transported  back  in  the  process  in 
such  a  manner  that  they  will  react  closer  to  the  surface.  Vandenker- 
ckhovo  thus  assumes  that  inside  the  fizz  zone  temperature  T^  is  exactly 
reached  at  the  limit  of  the  laminar  sublayer.  The  heat  transfer  equation 
therefore,  relates  +h°  sv'facc:  temperature,  and  thus  the  burning  rate, 
to  the  thickness  of  the  laminar  sublayer  which  depends  on  mean  flow 
velocity  through  fluid  dynamics  considerations. 

Like  Comer,"  he  approximated  the  flow  pattern  in  the  perforation 
by  the  well  known  theory  of  pipes,  in  which,  for  sufficiently  large 
Reynolds  numbers,  the  transition  from  laminar  to  turbulent  takes  place 
at  a  certain  distance  y  of  the  surface,  characterized  by  a  critical  value 
of  the  so-called  friction  distance  parameter  y  .  For  pipes,  this  critical 
value  ranges  from  5  to  10.  In  the  problem  under  consideration,  a  some¬ 
what  different  value  can  be  expected,  due  to  the  second  order  effect  of 
mass  addition  whose  velocity,  perpendicular  to  the  surface,  represents 
approximately  one  per  cent  of  the  mean  flow  velocity  in  the  region  of 
intere  st . 

This  infuence  car  be  taken  into  account  empirically  by  identifying 
the  threshold  velocity  U  which  is  the  mean  flow  velocity  below  which  no 
erosion  occurs,  as  the  precise  velocity  for  which  the  thickness  of  the 
fizz  zone  is  exactly  equal  to  that  of  the  laminar  sublayer. 
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In  the  numerical  example,  for  U  =  1 80  m/sec,  he  found  : 


2.85 


which  appears  reasonable. 


Above  the  threshold  velocity,’  turbulence  penetrates  the  fizz  zone, 
thus  reducing  its  effective  thermal  thickness  and  accordingly  increasing 
the  surface  temperature  and  the  burning  rate.  This  increase  can’be 
computed  easily  if  it  is  further  assumed  that  the  critical  friction 
distance  parameter  remains  constant.  > 


Despite  many  drastic  assumptions,  this  theory  has  been  used  to 
compute  a  value  of  the  erosive  constant  of  JPN  which  agrees  remarkably 
with  experimental  data.  Over  and  above  this  agreement,  this  work 
provides  a  plausible  physical  explanation  for  the  existence  of  the 
threshold  velocity  which  has  been  observed  in  several  experiences  with 
double  base  propellants  (see  figure  3  taken  from  Wirpress  [12]  but 
has  also  failed  to  appear  in  many  other- tests,  espeu  ally  with  composite 
propellants,  and  has  thus  been  considered  as  controversial  by  many 
authors  (an  attempt  to  explain  these  different  behaviours  is  made  later 
in  this  report,  in  section-  '4  —  B) . 


This  theory  still  suffers  from  several  drawbacks! 


-  the  influence  of  mass  addition  has  been  taken  into  account  in  a  most  - 
empirical  manner  while  the  influence  of  the  distance  from  the  fore  end 
of  the  grain  is  still  neglected. 


-  a  single  example  has  been  computed  and,  due  to  the  lack  of  experimental 
data,  no  attempt  has  been  made  for  studying  systematically  the  influence 
of  such  factors  as  burning  rate,  pressure  or  perforation  dimensions. 


-  finally  the  calculation  has  been  based  on  Rice  and  (Jinell  theoretical 
results  obtained  for  HES  401 6  double  base  propellant,  while  the  values 
of  the  threshold  velocity  and  erosive  coefficient  were  those  measured 
for  JPN. 
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Nevertheless  the  good  numerical  results  were  most  encouraging 
despite  the  great  numb?r  of  parameters  whose  values  were  only  approxim¬ 
ately  known . 

As  for  the  preceding  theories,  however,  this  agreement  with  ex¬ 
perimental  data  could  "be  fortuitous  or  could  result  from  parameter 
value  adjustment. 

The  main  purpose  of  the  present  report  is  to  check  the  validity 
of  Vandenkerckhove  theory  by  applying  it  to  composite  propellants  and 
especially  to  formulations  based  on  Ammonium  Perchlorate  which  have 
never  been  extensively  studied  before,  despite  some  preliminary  results 
for  Ammonium  Nitrate  propellants  [10,  11]'. 


3.  EXPERIMENTAL  DATA. 


A.  Burning  Rate. 


Figure  3»  taken  from  Wimpress,  is  typical  of  the  early  data  obtained 
for  particular  propellants  without  taking  apparently  into  account  such 
factors  as  pressure  or  motor  scale  [  12],  In  this  particular  case, -for 
J.P.N.,  the  existence  of  a  threshold  velocity  is  clearly  observed. 

At  our  best  knowledge,  the  first  systematic  experimental  study  of 
the  phenomenon  is  due  to  Green  who'  established  he  dependence,  between 
erosive  burning  and  strand  burning  rate,  for  both  double- -base 'and  compos¬ 
ite  propellants  [  1  3 ] . 

These  results  are  summarized  by  figure  4,  the  erosion  coefficient 
kj,  being  defined  by  the  relationship: 


r  =  rc  (  1  +  kM  G/G  ) 

*< 

r  being  the  burning  rate  without  erosion,  and  G/G  the  ratio  of  actual 
mass  velocity  PV  to  its  critical  value,  which  is  a  function  of  MACH 
Number  [ 1 ] . 

The  strand  burning  rate  r^  is  taken 'at  the  standardized  pressure 
of  70.3  kg/cm  although  the  actual  measurements  have  been  performed  at 
significantly  lower  pressure,  using  a  strongly  progressive  grain  (initial 
inner  radius  equal  to  1.5  cm.).  Green,  therefore,  has  not  taken  into 
account  the  influence  of  pressure  on  erosion  coefficient  as  well  as 
the  existence  of  the  threshold  velocity  which  was  unnoticed  in  most  of 
his  tests.  These  results  have  also  been  correlated  by  Geckler  with 
specific  heat  ratio  and  adiabatic  combustion  temperature  in  a  rather 
crude  manner  which  does  not  take  into  account  the  influence  of  ballistic 
modifiers,  oxidizer  granulation,  or  fuel-binder  nature  [  8j. 
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These  two  references  show  clearly  that  fast  burning  propellants 
are  le3S  sensitive  to  erosion  than  slow  burning  ones. 

The  best  investigation  of  erosive  burning  available  in  the  open 
literature  has  been  performed  in  Sweden  by  Lake  and  Karklund.*  [  141  . 
Their  results  are  represented  on  figure  5  and  discussed  in  details  in 
section  4  -  B.  It  is  seen  that  for  a  sufficiently  strong  erosion  (for 
G/G  >  0.2  )  the  severity  of  erosive  burning  increases  significantly 

with  pressure  as  shown  by  figure  6  which  gives  as  a  function 

of  pressure.  dG  /G 

Karklund  has  recently  obtained  new  and  still  unpublished  results 
with  both  pure  ammonium  perchlorate  and  composite  propellants  l  15  i  • 
These  tests  in  general  confirm  our  theoretical  analysis  (see  section  4). 


X  X  X  X 


An  increase  in  burning  rate  is  not  necessarily  always  associated 
with  erosive  burning;  in  fact  some  tests  suggest  that  for  small  erosive 
velocities,  below  50  m/sec,  the  burning  rate  of  certain  composite  for¬ 
mulations  may  be  slightly  decreased  below  r  [  13,  141  .  This  phenomenon 
is  still  unexplained;  it  may  perhaps  just  be  a  secondary  effect  of  the 
experimental  technique. 

On  the  other  hand,  Price  indicated  that  the  so-called  mesa  double 
base  propellants  usually  exhibit  decreases  in  burning  rate  during  oscillat¬ 
ory  burning  l  16]  .  This  observation  suggests  that  similar  decreases  may 
result  from  erosive  bruning. 


B.  Influence  of  Grain  Design  and  Kouor  Scale. 


Experimental  evidence  shows  that  erosive  burning  is  also  influenced 
by  grain  design  and  possibly  by  port  dimension. 

In  an  attempt  to  correlate  the  erosive  coefficient  k:T  with  grain 
design,  Fenech  and  Billheimer  introduced  a  configuration  factor  y  which 
is  a  measure  of  grain  complexity: 


A  4  “  A 

P  P 

A  being  the  flow  area  of  a  circular  perforation  having  the  same  perimeter 
Scthan  the  actual  port  whose  area  is  A^  f  1 7  1  - 


It  is  believed  that  a  similar  investigation  is  described  by  Saderholn  in 
a  Thiokol  special  report,  "Erosive  Burning  Study  of  TP  -  H  8041  Propellant" 
(Confidential) . 
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Using  Green's  relationship,  these  authors  computed  the  values  of  k^r 
which  permit  to  fit  best  the  pressure-time  traces  obtained  by  firing  the 
eight  grains  represented  by  figure  7.  The  results  of  this  computation  and 
the  comparison-  with  the  data  obtained  by  Green  are  represented  on  figure 
8  from  which  it  is  seen  that  in  most  cases  the  erosive  coefficients  are 
significantly  smaller  than  those  measured  by  Green  although  in  two  cases 
they  fall  almost  exactly  on  his  average  line. 

These  erosion  coefficients  k^  are  also  represented  as  a  function  ■ 
of  the  configuration  factor  y  on  figure  9  which  suggests  that  erosive 
burning  increases  with  grain  complexity,  the  least  squares  line  equation 
being* 

0.38. 

k.,  =  0.21  V 

M 

i 

The  eight  tested  grains  were  differing  much  in  configuration 
factor  y  ,  in  port  to  throat  area  ratio  H,  and  in  standard  burning 
rate  rQ  std,  at  70.3  kg/cm  ,  as  indicated  by  figure  7. 

It  is,  however,  most  difficult  to  draw  definite  conclusions  from 
these  results'  since  the  mean  chamber  pressure  and  the  port  area 
(representing  the  scale)  of  these  motors  have  not  been  reported.  It  is 
most  probable  that  these  two  factors  have  a  significant  influence  on 
k^  and  under  the  assumption  of  similar  chamber  pressures  figure  8 
suggests  that  the  experiments  of  Fenech  and  Billheimer  have  been  performed 
with  full  scale  motors  of  much  larger  dimensions  than  the  small  test 
rockets  used  by  Green. 

It  must  be  emphasized,  however,  that  the  interpretation  of  these 
results  is  most  delicate  in  terms  of  the  erosion  coefficient  k^.  defined 
by  an  empirical  equation  which  does  not  take  into  account  the  possible 
existence  of  a  velocity  threshold.  For  instance,  it  is  not  surprizing 
that  grain  B  of  figure  7,  whose  port  throat  area  ratio  is  equal  to  10, 
is  characterized  by  the  lowest  value  of  k^.  Indeed,  under  such 
conditions  of  rather  low  erosion,  it  is  probably  operating  near  the 
threshold  velocity  and  Green's  equation  can  hardly  be  expected  to 


>(  =  11.0  '  H  =2.25  >(=1.0  H  =  10.0 

r  =  0.790  cm/sec  r  .  .  =  0.739 
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represent  accurately  the  actual  phenomenon. 

The  results  obtained  by  Fenech  and  Billheimer  are  average  values 
computed  in  the  assumption  of  an  equal  burning  rate  along  the  perimeter 
of  each  cross-sectional  area.  In  fact  Dickinson,  Jackson  and  Odgers 
reported  higher  burning  rates  on  the  tips  of  the  arms  of  certain  grain 
design,  than  on  the  rounded  portions  and  sides  of  the  arms  l  1 8  1  . 

These  authors  have  shown  that  this  difference  in  burning  rates  decreases 
when  the  grain  design  is  such  that  the  local  ratio  of  burning  perimeter 
to  available  port  area  for  gases  does  not  vary  much  from  one  point  to 
another,  in  a  cross-section.  They  attributed  possibly  the  phenomenon 
to  cross-flow.  Another  explanation  can  be  found  in  the  influence  of 
local  port  scale  which  can  be  assumed  to  cover  both  the  influence  of 
motor  scale  (overall  effect)  and  that  of  grain  design  (local  effect). 

A  quantitative  theory  in  this  direction  is  given  in  section  5. 


C.  Experimental  Techniques. 


Without  entering  into  the  detailed  description  of  the  experimental 
techniques,  it  is  essential  to  emphasize  that  measurements  on  ero&ive 
burning  can  be  obtained,  by  two  ways. 

In  the  first  technique  the  burning  rate  measurement  is  made  on 
the  main  grain,  either  by  interrupting  the  combustion  [12,  13]  through 
suitable  probes  embedded  in  the  propellant.  [181  or  by  analysis  of  the 
engine  performance  curve  [171  ♦  Figure  10  represents  the  experimental 
motor  used  by  Green  [  13 1. 

In  the  second  technique,  the  test  sample  is  separate  from  the 
main  charge,  and  in  general  located  in  a  relatively  long  pipe  leading 
to  the  nozzle  [  14,  1.51  •  In  some  cases,  however,  the  nozzle  itself,  or 
part  of  it,  is  made  of  propellant  (5,  6,  7,  141-  Figure  11  represents 
the  different  arrangements  used  by  Lake  and  Marklund  which’ have  n_ 
reported  a  significant  difference  between  the  measurements  made  with 
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tablets  and  strips  C  1 4- »  15  3  • 

It  i3  probable  that  the  selected  experimental  technique  somewhat 
influences  the  measurements.  Indeed  in  the  first  method, as  pointed 
out  before,  the  length  to  diameter  ratio  of  the  perforation  is  seldom 
sufficiently  large  for  neglecting  the  two-dimensional  nature  of  the 
flow  field  which  i3  influenced  by  the  distance  from  the  fore  end  of 
the  grain  and  by  mass-addition. 

On  the  other  hand,  in  the  second  method  using  a  sufficiently  long 
pipe,  a  fully  developed  turbulence  can  bo  assumed,  nevertheless  the 
influence  of  the  interaction  between  main  flow  and  mass-addition  due 
to  the  test  sample  can  hardljr  be  assessed. 

In  fact  the  flow  fields  obtained  by  the  two  methods  are  somewhat 
different  and  it  is  not  surprizing  that  the  results  obtained  by  the 
second  technique  yield  in  many  instances  somewhat  larger  values  of 
erosive  burning  than  those  measured  on  the  main  charge  (see  section  3  -  B) . 
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4.  THEORETICAL  STUDY  OP  EROSIVE 


BURNING  OP  AMMONIUM 


PERCHLORATE  PROPELLANTS 


A.  Burning  Mechanism  of  Ammonium  Perchlorate  Propellants. 


Before  proposing  a  model  of  erosive  "burning  of  ammonium  perchlorate 
propellants,  it  is  necessary  to  study  in  some  details  the  "burning 
mechanism  of  these  propellants. 

The  problem  is  presently  rather  well  understood,  owing  to  the 
recent  results  obtained  by  professor  M.  Summerfield  and  his  group  at 
Princeton  University  [  1 9 1  who  have  confirmed  the  theoretical  predictions 
previously  proposed  by  Jaumotte  and  Vandenkerckove  l  2.0l  from  which  the 
following  analysis  is  partially  drawn. 

Much  information  can  be  obtained  simply  by  comparing  the  burning 
rate  of  propellants  to  that  of  pure  ammonium  perchlorate,  as  reported 
by  Friedman  and  his  coauthors  t  21 ,  22  ]  .  The  burning  rates  of  four 
propellants  using  polystyrene  as  fuel-binder  and  differing  only  by 
mixture  ratio  and  granulation  are  taken  from  Summerfield  and  Marklund’s 
works  f  14;  23  ]  .  The  comparison  given  by  figure  12,  shows  that  the 
burning  rate  of  the  propellant  can  be  higher  or  lower  than  that  of 
pure  oxidizer,  depending  on  pressure,  mixture  ratio  and  granulation; 
it  can  be  summarized  as  follows: 

a  -  fojr-coarse  granulations  and  low  oxidizer  contents,  the  propellant 
burning  rate  is  higher  at  low  pressure  and  smaller  at  high  pressure 
than  that  of  pure  perchlorate . 

b  -  for  fine  granulations  and  high  oxidizer  contents,  the  burning  _ute 
is  higher  at  all  pressures  in  the  considered  range-, 
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c-  below  approximately  50  kg/cm  ,  the  combustion  index  n  of  the  propellant 

(as  given  by  the  classical  law  r  =  apn  )  is  definitely  smaller  than 

that  of  pure  perchlorate  which,  apparently, „ becomes  higher  than  unity 
2 

below  40  kg/cm  . 

d-  For  pure  perchlorate,  Friedman  indicates  a  limit  of  flammability  at 
2 

about  22  kg/cm  while  Summerf ield .has  been  able  to  measure  steady 
propellant  combustion  down  to  one  atmosphere.  This  fact  suggests  that  the 
burning  mechanisms  are  somewhat  different  in  the  low  pressure  range, 
e-  Finally  the  near  constant  combustion  index  n  exhibited  by  most  propellants 
in  the  full  pressure  range  does  not  suggest  that  the  phenomenon  which 
controls  the  propellant  burning  rate,  namely  some  kind  of  interaction 
between  the  oxidizer  and  the  fuel,  is  basically  different  above  and 
below  the  burning  rate  of  pure  perchlorate. 

X  X  X  X 

It  seems  reasonable  to  assume  that  in  the  propellant  a  single  oxidizer 
particle  burns  at  least  as  fast  as  pure  perchlorate  at  the  same  pressure. 
Indeed  the  decomposition  zone  of  the  oxidizer,  being  of  the  order  of  one 
micron  or  less,  is  generally  much  thinner  than  the  particle  dimensions 
and  the  effect  of  the  interface  can  thus  be  neglected.  The  express:  on 
"at  least  as  fast"  has  purposedly  been  used  since  the  burning  rate  of 
pure  perchlorate  appears  as  a  lower  limit  and  heat  can  conceivably  be 
brought  from  the  main  flame  zone  back  to  the  oxidizer  for  increasing  its 
temperature  and  deflagration  rale. 

It  appears  therefore ,  that  in  the  region  situated  below  the  deflagration 
curve  of  pure  ammonium  perchlorate,  the  burning  rate  of’ the  oxidizer  is 
faster  than  the  average  burning  rate  and  that  fuel  peaks  are  most  probably 
protruding  the  surface . 

On  the  other  hand,  above  the  deflagration  curve  of  pure  ammonium 
perchlorate,  no  clear  conclusions  can  be  drawn  at  this  point  since  the 
lower  limit  of  the  oxidizer  deflagration  rate  in  the  propellant  is 
only  known. 
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If  the  decomposition  flame  is  the  only  process  of  importance,  the 
perchlorate  particles  will  protrude  above  the  surface  of  the  faster  pyroly- 
zing  fuel;  if  significant  heat  transfer  exists  between  the  diffusion 
flame  and  the  surface,  however,  both  situations  can  prevail  and  we  must 
rely  on  experiments. 

In  fact,  early  experiments  have  shown  that  above  the  deflagration 
curve  of  the  pure  oxidizer,  at  atmospheric  pressure,  large  ammonium 
perchlorate  crystals  are  protruding  from  the  main  surface  [3,  231. 

Later  experiments  provided  much  more  precise  information.  Using 
high  speed  cinematography,  Barr&re  observed  a  decrease  in  oxidizer 
particle  protrusion  with  increasing  pressure  [  24].  The  interpretation 
of  these  excellent  cinematographic  records,  for  determining  the  surface 
structure,  remains  delicate,  especially  in  the  high  pressure  range. 

After  similar  attempts,  Summe rfie Id' s  group,  therefore,  developed 
a  method  whereby  small  pieces  of  burning  propellant  can  be  extinguished 
suddenly  by  rapid  evacuation  of  the  burning  vessel  due  to  burst  of  a 
rupture  disk  l 181.  This  method  is  believed  to  leave  the  propellant 
surface  intact  and  unaffected  by  the  rapid  pressure  drop,  for  microscopic 
study  of  the  structure.  The  pictures  taken  by  this  method  show  that  at 
low  pressure  the  oxidizer  crystals  project  above  the  fuel  surface.  As 
the  pressure  increases,  the  crystals  become  flush  and  eventually  recede 
deeply  into  the  fuel,  as  schematized  by  figure  13.  Very  finely  ground 
oxidizer  particles,  however,  may  be  thrown  away  from  the  surface  as  soon 
as  they  are  reached  by  the  flame . 

For  our  present  purpose  several  conclusions  can  be  drawn  from  these 
observations . 

a  -  In  the  low  pressure  range,  below  30  kg/cm  approximately,  the  burning 
rate  of  the  propellant  is  much  higher  than  that  of  the  oxidizer,  and 
most  probably  some  heat  must  be  transferred  from  the  main  reactions 
between  oxidizer  and  fuel  vapors,  back  to  the  surface, 
b  -  It  seems  doubtful,  however,  that  this  heat  transfer  is  homogeneous 
all 


over  the  surface . 
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Indeed,  at  the  considered  pressure,  according  to  Summerfield,  the  char¬ 
acteristic  chemical  time  of  the  reaction  between  oxidizer  and  fuel  va¬ 
pors  is  much  shorter  or  at  least  of  the  same  order  of  magnitude  than  the 
characteristic  diffusion  time  £18]  .  Highly  exothermic  reactions,  there¬ 
fore,  most  probably  take  place  very  close  to  the  surface  at  the  interface 
of  the  oxidizing  and  reducing  streams.  In  this  region,  only  a  small 
fraction  of  the  gas  is  involved  hut  this  may  he  sufficient  for  modifying 
the  surface  temperature  locally  at  the  interface  which  has  its  ovn 
regression  rate. 

It  appears  therefore,  that  the  regression  rate  of  the  interface 
between  oxidizer  and  fuel  is  of  particular  importance  in  controlling  the 
burning  rate  of  the  propellant,  in  the  low  pressure  range.  This  conclusion 
is  consistent  with  the  effect  of  fuel  nature,  and,  overall,  of  cixturo 
ratio  and  granulation.  Indeed  a  larger  oxidizer  percentage  and  a  finer 
granulation  "both  increase  the  length  of  the  rate  controlling  interface 
and  the  number  of  crystals  reached  per  unit  time,  thus  reducing  the 
distance  actually  travelled  by  the  interface  regression  per  unit  length 
burned  perpendicular  to  the  average  surface. 

c  -  It  must  be  emphasized,  however,  that  in  the  whole  pressure  range  the 
regression  rates  of  the  fuel,  oxidizer  and  interface  are  not  much  diffe¬ 
rent.  Ir.  particular  in  the  low  pressure  range  the  interface  to  oxidizer 
regression  rate  ratio  cannot-  be  much  higher  than  1  .57  unless  ammonium 
perchlorate  particles  are  rooted  up  as  represented  very  schematically  hy 
figure  14  drawn  in  the  oversimplified  assumption  that  the  regression  of 
a  spherical  particle  is  unaffected  by  the  interface. 

In  fact  the  actual  phenomenon  is  almost  certainly  much  morb  complex} 
the  interface  regression  rate  may  be  influenced  hy  other  factors  such  a3 
particle  size  and  time,  hut  the  general  conclusion  of  similar  regression 
rates  still  remains  valid. 

d  -  In  the  high  pressure  range,  above  approximately  60  kg/cm  ,  the 
burning  rate  of  all  the  propellants  whose  tests  have  been  reported  in 
the  open  lit:rature  is  never  much  higher  than  the  regression  rate  of 
the  oxidizer,  at  least  in  the  absence  of  ballistic  modifiers.  It  can, 
however,  be  much  lower  but  it  is  in  fact  often  rather  similar  l  23]  • 


umu 
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This  suggests  that  in  the  high  pressure  range  the  main  diffusion  flame 
does  not  play  anymore  a  significant  role  in  transfering  heat  towards  the 
surface  since  otherwise  the  propellant  "burning  rate  would  "be  consistently 
larger. 

It  is  not  our  purpose  to  discuss  much  further  the  complex  "burning 
mechanism  of  ammonium  perchlorate  propellants;  the  preceding  considerations, 
however,  are  necessary  for  studying  erosive  "burning  of  such  propellants. 


B.  Erosive  Burning  of  Ammonium  Perchlorate  Propellants. 


Our  work  relies  heavily  on  the  data  obtained  by  Lake  and  Marklund, 
and  represented  on  figure  5  £  14 1.  These  results  exhibit  a  most  interest¬ 
ing  behaviour. 

At  high  pressure  (100  kg/cm  )  the  burning  rate  without  eroBion  of 
the  three  considered  propellants  depends  in  the  normal  way  on  mixture 
ratio  and  fuel  nature,  the  main  particle  diameter  remaining  constant 
(24  -  30  0  in  the  experiments. 

Surprizingly  enough,  however,  in  the  presence  of  a  sufficiently 
strong  etosion,  above  100  m/sec,  the  influence  of  mixture  ratio  and 
fuel  nature  quickly  disappears,  and  all  three  propellants  exhibit  an 
identical  burning  rate  -  versus  -  gas  velocity  relationship. 

It  appears  therefore,  that  the  burning  rate  without  erosion  can 
be  viewed  as  an  abnormal  behavior  and  that  in.  the  presence  of  erosion,  a 
classical  burning  lav  of  the  form  : 

r  =  r^  multiplied  by  a  function  of  erosion  severity 
is  misleading  since  the  reference  rate  rQ  does  not  influence  appreciably 
the  burning  rata-versus-velocity  relationship. 

On  the  other  hand,  at  medium  and  low  pressures  (50  and  20  kg/cm  ) 
this  behaviour  does  not  hold  true  anymore. 


In  these  experiments,  indeed,  the  influence  of  fuel  nature  remains 
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in  the  presence  of  erosion,  while  the  influence  of  mixture  ratio  still 

2 

disappears  at  50  kg/cm  and  insufficient  data  leaves  the  question  open 
2 

at  20  kg/cm  . 


X  X  X  X 

In  an  attempt  to  explain  these  phenomena,  the  burning  rates  without 
erosion  of  the  three  propellants  used  by  Marklund  and  Lake  have  also  been 
plotted  on  figure  12.  The  correlation  between  figures  5  and  12  is 
interesting.  At  high  pressure,  when  the  strand  burning  rates  of  the  low 
mixture,  ratio  propellants  tested  in  Sweden  are  below  the  deflagration 
curve  of  pure  ammonium  perchlorate,  the  influence  of  composition  disappears 
with  erosion  while  at  medium  and  low  pressures,  when  the  burning  rate 
of  the  faster  composition  is  above,  the  influence  of  composition  subsits, 
at  least  partially.  The  fact  that  at  50  kg/cm  the  influence  of  fuel 
nature  subsists  in  the  presence  of  erosion  while  the  influence  of  mixture 
ratio  still  disappears  suggests  that  a  distinction  must  be  made  between  the 
effects  of  fuel  nature  and  mixture  ratio  on  burning  mechanism. 

X  XXX 

Furthermore  one  can  expect  that  the  controlling’ factor  of  these 

» 

phenomena  is  not  the  relative  magnitude  of  the  propellant  and  oxidizer 
regression  rates,  but  rather  the  surface  structure. 

Taking  for  granted  that  erosive  burning  results  primarily  from  the 
increasing  heat  transfer  which  takes  place  just  above  tHe  surface,  due 
to  the  flow,  the  following  physical  model  can  be  proposed.  With  the 
high  rugosity  characterising  the  burning  surface  one  can  hardly  assume 
the  existence  of  a  normal  laminar  sublayer  above  the  surface  protrusions, 
but  it  is  reasonable  to  assume  that  these  peaks  are  eroded  much  more 
severely  than  the  rest  of  the  surface  which  then  becomes  smoother  and 
smoother  with  increasing  erosion. 

a  -  At  low  pressure  when  these  protrusions  consist  in  oxidizer  particles, 
their  erosion  does  not  influence  appreciably  the  burning  rate  .intil 
the  surface  becomes  sufficiently  flat  and  an  apparent  "threshold 
velocity"  exists  below  which  no  significant  erosion  takes  place. 
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In  the  presence  of  a  strong  erosion  the  fuel  nature  still  influences 
the  turning  rate,  a  phenomenon  which  confirms  the  conclusion  that  the 
regression  of  the  oxidizer  particles  is  markedly  accelerated  by  heat  transfer 
from  at  least  part  of  the  diffusion  flame. 

b  -  At  medium  pressure  some  fuel  peaks  have  probably  appeared, th$>  threshold 
velocity  has  almost  disappeared  except  for  the  higher  mixture  ratios 
but  the  fuel  nature  still  influences'  somewhat  the  burning  rate  under 
severe  erosion  conditions.  • 

c  -  At  high  pressure,  the  fuel  peaks  which  protrude  from  the  surface  are 
first  eroded.  Due  to  the  high  oxidizer  content  their  erosion  results 
in  new  oxidizer  crystals  being  reached  by  the  flame  sooner  than  without 
erosion. 

In  this  case,  no  threshold  velocity  exists  and  figure  5  shows  that 
the  steepest  increase  in  burning  rate  with  gas  velocity  corresponds  to 
the  largest  difference  between  erosionless  burning  rate  and  deflagration 
rate  of  pure  perchlbrate  (largest  fuel  protrusions). 

Furthermore  it  has  been  concluded  previously  that  at  high  pressure 
no  significant  heat  from  the  main  diffusion  flame  reaches  the  oxidizer 
particles.  These  conclusions  are  confirmed  by  the  experimental  fact  that 
the  influence  of  fuel  nature,  mixture  ratio  and  probably  granulation 
disappears  when  erosion  becomes  sufficiently  severe,  thus  indicating  that 
the  decomposition  of  th«  oxidizer  crystals  becomes  the  single  'controling 
factor  of  the  burning  rate,  the  surface  being  assumed  to  be  sufficiently 
smooth  (without  erosion  the  burning  rate  is  also  controlled  to  some  extent 
by  granulation  and  mixture  ratio,  a  fact  which  could  perhaps  be  explained 
by  the  probability,  for  the  flame  to  pass  from  one  exidizer  crystal  to 
another) . 


Fhen  the  surface  is 'flattened  by  erosion,  the  oxydizer,  fuel  and 
interface  regression  rates  become  equal.  At  high  pressure  the  problem 
can,  therefore,  be  treated  as  if  the  propellant  is  replaced  by  pure  ammonium 
perchlorate  (possibly  catalyzed  by  suitable  ballistic  modifiers). 


'A 


X 


X 
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For  checking  the  validity  of  this  physical  model,  it  is  possible 
to  compute  the  slope  of  the  experimental  ~  versus  -  rQ  ourve  obtained 
by  Green  (see  figure  4  taken  from  f 1 3 1  )  which  proposed  to  represent 
erosive  burning  by  the  empirical  expression  t 

r  ■  r  (1  e  k,.  G/G*) 

0 

rQ  being  the  burning  rate  without  erosion,  k^  the  erosive  coefficient 
and  G/G*  the  ratio  of  aotual  mass  velocity  f>  V  to  its  critical  valued]. 

The  following  computation  is  carried  out  at  the  reference  pressure 
of  70  kg/cm2  (about  1.000  psi)  which  can  still  be  considered  as  a  rel¬ 
atively  high  pressure. 


At  that  particular  pressure,  interpolation  of  the  data  obtained  by 
Lake  and  Marklund  can  be  represented  approximately  by  a  straight  line 
r  -  versus  —  G/G*  when  rQ  is  equal  to  0.65  om/sec,  and  r  ■  1.97  on/soc 
for  o/G*  =  1,  as  schematized  on  figure  15.  (see  reference  [1  ]  and  its 
figure  7  for  complete  treatment  of  this  problem  of  internal  ballistics, 
taking  into  ao count  that  .j 


(  G 


h  * >/  -0.529  for  y  e  1*24.) 


The  empirical  equation,  therefore,  givasi 


1  .97 

0.63“ 


-  2.035 


This  value  is  much  higher  than  the  one  measured  by  Green,  for 
that  particular  burning  rate,  i.e. 


k y.  *  0.98 

and  the  two  series  of  experiments  can  only  be  correlated  through  the  use 
of  a  rather  large  coefficient  » 

2 ...  0  >  5  /  0.98  -  2,080 

It  is  "tnl  level  'hat  This  discrepancy  can  be  explained  as  folxovo  t 


BURNING  RATE  (cm/sec) 


FIG. 15 
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1 /  The  experimental  techniques  are  different  as  emphasized  in  section  3. 

2/  Although  Green  correlated  his  results  "by  using  the  reference  burning 

rates  at  'JO  kg/cm  ,  his  measurements  have  actually  "been  made  at  a  much 
smaller  pies  sure,  at  least  at  the  "beginning  of  the  "burning  time  when  eros¬ 
ive  burning. is  the  largest  (  a  strongly  progressive  grain  design  was 
used.) . 

Under  these  conditions,  an  appreciably  smaller  erosive  coefficient 
is  consistent  with  the  data  obtained  by  Lake  and  Marklund. 

In  addition  to  a  correlation  coefficient  equal  to  2.080,  if  it  is 
assumed  that  Green's  measurements  have  been  obtained  near  G/G*  «  1, 
the  slope  of  the  -  versus  -  r"  curve  can  be  easily  computed  in  the 
neighbourhood  of  the  chosen  point  of  reference  (r  »  0.65  cm/sec  ,  k^  =  O.98). 

Indeed  the  physical  model  predicts  that  under  a  severe  erosion,  the 
burning  rate  is  independent  of  r  and  equal  to  1  .97  cm/sec  in  the  case  Tinder 
consideration. 

The  empirical  equation,  therefore,  gives  for  q  =  1  cm/sec  t 

1  .97  -  1  (  1  +  1  .080  kjj)  . 

or  k^  =  O.466 

Similarly,  for  rQ  =  0.5  and  0.2  cm/sec,  one-  gets  respectively 
=  1.415  and  4.250. 

These  results  are  compared  with  Green's  data  on  figure  16.  The 
correlation  is  fairly  good,  especially  if  the  dispersion  of  the  experimental 
points  is  taken  into  account.  . , 

Actually  Green's  data  have  been  obtained  for 
G/O*  ^  1 

and  a  slightly  higher  correlation  coefficient  would  have  to  be  used 
(2.12  for  G/G*)  resulting  in  a  somewhat  larger  slope  of  the  computed  curve. 


39 


It  must  be  emphasized,  however,  that  the  equation  which  defines 
k^  is  only  an  approximate  and  empirical  expression  which  "barely  represents 
the  actual  phenomenon  in  the  low  velocity  range  (see  the  difference 
between  the  plain  and  dotted  lines  on  figure  15)  and  that  it  is,  consequent¬ 
ly,  difficult  to  carry  the  comparison  .  beyond  this  point. 

Furthermore  it  must  be  kept  in  mind  that  Green's  measurements  were 
taken  at  relatively  moderate  pressure,  in  a  range  in  which  the  fuel  nature 
and  possibly  other  factors  are  still  of  some  importance  under  severe 
erosion. 

In  view  of  the  drastic  assumptions  made,  the  preceding  correlation 
does  not  prove  rigorously  the  validity  of  the  proposed  model  $  it  is 
nevertheless  encouraging. 
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C.  Mathematical  Model  of  Erosive  Burning. 


£  Solid-Phase  Decomposition. 

The  following  method"-  is  similar  to  the  one  proposed  in  reference  [  9  I 
for  computing  the  erosive  Burning  of  double • base  propellants,  the  so- 
called  fizz  zone  being  replaced  by  the  oxidizer  decomposition  zone  at 
the  end  of  which  temperature  T.  is  obtained.  As  emphasized  previously 
this  approach  is  valid  only  at  relatively  high  pressure,  above  5®  kg/cm  , 
in  a  region  where  no  significant  heat  transfer  from  the  main  flame  zone 
takes  place  towards  the  surface,  thus  permitting  to  compute  temperature 
by  assuming  an  adiabatic  decomposition  reaction. 

Under  severe  erosion  ,  therefore,  the  propellant  can  be  represented 

by  a  fiat  surface  of  pure  ammonium  perchlorate  and  the  solid-phase 

decomposition,  which  is  considered  by  the  majority  of  the  authors  as  the 

rate  controlling  step  of  the  combustion,  is  given  in  terms  of  the  surface 

temperature  T  by  the  Arrhenius  type  relationship  |  3,  8,  25  1  : 
s 

r  «=  B  exp  (  -  E  /  Ro  Ts  )  (l) 

with  r  the  hurtling  rate 

B  the  preexponentiu.1  fav/toj. 

E  the  activation  energy 

and  Rq  the  universal  gas  constant  ( 1  . 9 85  cal/g.mole.  °K) 

Certain  authors ^however,  use  a  slightly  different  expression!  26,  27! « 

r  =  B'  Tg  exp  (  -  E  /  Ro  Ts  )  0') 

which  in  practice  gives  similar  results. 

f>  Heat  Transfer  across  the  Decomposition  Zone. 


It  is  assumed,  hereunder,  that  the  heat  transfer  across  the  '-as-phase 
decomposition  zone  is  due  primarily  to  gas  conductivity  whereas  the 
influence  of  turbulence  and  radiation  can  be  neglected  to  a  first  approxim¬ 
ation. 
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At  a  distance  y^  away  from  the  burning  surface,  the  chemical  reactions 
characterizing  the  decomposition  zone  eve  completed,  and  in  the  absence  of 
any  heat  transfer  from  the  flame  zone  temperature  is  obtained. 

In  order  to  obtain,  as  a  function  of  y,  the  temperature  profile  within 
the  decomposition  zone,  the  degree  of  completion  of  the  reaction  at  each 
distance  y  from  the  surface  is  needed.  This  could  be  done  through  chem¬ 
ical  kinetics  considerations  along  the  lines  used  by  NACHBAR  and  others  t  25]. 
This,  however,  would  increase  very  much  the  complexity  of  the  problem  and, 
overall,  introduce  many  parameters  whose  numerical  values  are  not  accurately 
known,  thus  decreasing  the  confidence  level  in  the  results. 

In  accordance  with  several  authors  (Rice  and  Ginell,  Parr  and  Craw¬ 
ford,  Geckler),  therefore,  it  is  assumed,  hereunder,  that  to  a  first 
approximation  the  total  heat  release  takes  place  at  once  at  the  end  of  the 
decomposition  zone,  in  the  plane  y  =  y^  l  8  ... 

In  relative  motion,  the  oxidizer  being  fixed  at  the  burning  surface, 
no  heat  source  is  therefore  encountered  for 


0  <  y  <  y, 

and  in  steady  state  the  conduction  equation  can  bo  written: 


k 


dt 

dy 


Ppr 


c  (  T  -  T'  ) 
P  s' 


with  k  the  thermal  conductivity, 

pp  the  propellant  specific  mass, 

P  r  the  flow  rate  per  unit  area  of  the  surface, 
P 

and  c  the  specific  heat  at  constant  pressure. 


(2) 


This  realtionship  expresses  the  equality  of  the  heat  flux  by  conduc¬ 
tion  towards  the  surface  and  of  the  heat  flux  leaving  the  surface  '  o  the 

form  of  enthalpy.  The  term  c  T'  represents  the  energy  level  of  reference 

P  s 

whose  expression  is  established  below. 
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By  introducing  mean  values  of  the  specific  heat  c  and  thermal 
_  P 

conductivity  k,  the  preceding  equation  can  he  readily  integrated  between 

the  burning  surface  (y»0,  T-T)  and  the  end  of  the  decomposition 

zone  (  y  =  y1 ,  T  =  T1 ) .  It  yields  t 


T  -  T' 
s  s 

m  _  T, 
1  S 


-  P  r  — -r— ■  7i 


The  value  of  the  energy  level  of  reference  c^  T^  can  be  computed 
as  follows:  on  the  burning  surface,  for  y  =  0,  a  gaseous  phase  is  produced 
at  temperature  T  ,  and  the  passage,  at  that  particular  temperature,  from 
solid  to  gas,  can  be  represented  by  a  reaction  of  the  form: 


Z  n.  s. 
l  i  i 


Z  m  .  C-  . 
1  1  1 


S.  and  G.  representing  respectively  the  solid  and  gaseous  species;  n.  and 
m.  the  corresponding  numbers  of  moles.  The  quantity  of  heat  Q  sinvolved 

1  'S 

in  this  generally  endothermic  process  is  x 


T 

2m.  AH 

ili 


f  ni  \  A  Hf 


|  i  H ^  J  representing  the  heat  of  formation,  at  temperature  Tg,  of 
the  considered  species. 


The  term  c  T'  has  been  introduced  as  the  energy  level  of  reference 
P  3 

per  unit  mass.  The  heat  reaching  the  burning  surface  serves  for  heating 
the  solid  propellant  from  its  initial  temperature  T^  up  to  T^  and  then 
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for  achieving  the  phase  change,  a  process  which  is  generally  accompanied 
by  eome  chemical  reactions,  as  indicated  by  equation  (4).  The  energy 
balanco  per  unit  mass  can,  thus,  he  written  t 


C  (  T8  -  Tt  )  ♦  «, 


T'  ) 

s 


—  being  the  mean  specific  heat  of  the  solid  phase  and  the  specific  heat  cp 
being  assumed  to  remain  constant  in  the  decomposition  zone. 

The  heat  of  formation  is  not  readily  tabulated  at  temperature  Tg  hut 
rather  at  the  standard  temperature  Tq  =  25  °C,  yielding  » 


%*  -  -  0S  (  T8  -  T0  )  +  V  (  )  (7 

Combining  equations  (6)  and  (7)  gives  : 

r  (t„  -  t;  )  -  ^  <To  -  h  >  (a 

from  which  the  value  of  can  be  obtained. 

The  term  (T  -  T'  )  represents  the  fictitious  temperature  drop  which 
would  be  necessary  Jor  vaporising  aaiabatically  the  solid,  assuming  a 
cons  taut  mean  value  of  the  specific  heat  at  constant  pressure. 


For  an  endothermic  surface  reaction  ,  this  term  is  positive  and  it  is 
thus  normal,  although  somewhat  surprising,  that  T;.can  take  large  negative 
values  since  it  is  not  an  actual  temperature  but  rather  represents  a 
quantity  of  heat  with  its  sign. 
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In  the  presence  of  erosion  due  to  a  flow  parallel  to  the  propellant 
surface,  turbulence  takes  place  in  the  combustion  layer,  increasing  the 
temperature  level  and  decreasing  the  thickness  of  the  decomposition  zone. 
The  surface  temperature,  therefore,  increases  together  with  burning  rate, 
and  equation  (l')  can  be  written: 


T 

r  /S\P  E/1  1  \  /  \ 

—  =  (ip —  )  exP  “  ^  T -  ”  T —  '  (9) 

0  so  0  so  S 


subscript  indicating  the  absence  of  erosion. 

Within  the  decomposition  zone,  temperature  T^  is  obtained  at  a  point: 

y  '  <  y  > 

1 

where  the  decomposition  reactions  are  possibly  not  yet  entirely  completed. 
The  conduction  equation  can  be  written  directly  between  that  point  y.j  and 
the  burning  surface,  across  a  region  in  which  the  flow  can  be  assumed  to 
remain  laminar. 


Equations  (9)  and  (10)  permit  to  compute  T  and  r/r  corresponding  to 

s  o 

each  value  of  the  effective  decomposition  zone  thickness  : 

-  >r1 

v  -  Fluid  Dynamics. 


The  dependence  of  erosive  burning  on  the  effective  decomposition  '  ~ne 


thickness  y’  ,  at  the  end  of  vrhich  temperature  T^  is  obtained, 
sxipplemenied  by  a.  relationship  between  y  ^  and  the  mean  flow  v 
(or  Kach  number,  or  mass  vel...-3  y;  parallel  to  the  propellant 


must  be 

1  +TT  v 

-j 

surface  . 


i 
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Within  the  flame,  the  picture  of  separate  zones  is  probably  fairly  good 
when  it  is  applied  to  the  case  of  an  end-burning  strand  where  the  gas  flow 
is  approximately  laminar  (some  turbulence,  however,  appears  to  take  place 
at  the  interface  between  fuel  and  oxidizer  vapors)  and  one-dimensional, 
directed  normally  away  from  the  surface. 

In  the  much  more  complex  case  of  a  significant  velocity  parallel  to  the 
surface,  one  can  raise  the  question  of  the  validity  of  this  model. 

It  has  been  assumed,  hereunder,  that  the  presence  of  a  strong  velocity 
gradient  in  the  boundary  layer  does  not  alter  appreciably  this  picture,  in 
the  laminar  sublayer  where  all  transport  of  heat  and  chemical  species  is 
accomplished  by  movements  on  a  molecular  scale. 

In  the  turbulent  region,  however,  eddies  play  an  important  role  in  the 
transport  phenomenon  and  it  is  clear  that  the  picture  of  separate  zones 
does  not  hold  anymore.  Nevertheless  one  can  expect  that  this  flame  structure 
still  holds  in  the  transition  sublayer  where  turbulence  is  not  yet  fully 
developed  and  where  the  dimensions  of  the  addies  are  much  smaller  then  the 
thickness  of  a  zone. 

Furthermore  one  can  expect,  as  a  working  assumption,  that  inside  the 
decomposition  zone  temperature  T.  is  obtained  exactly  at  the  limit  of  the 

I 

laminar  sublayer.  Indeed  turbulence  in  the  upper  layers  of  the  flame  is 
not  yet  fully  developed,  its  influence  is  not  only  to  increase  the  thermal 
conductivity  hy  transporting  hot  gases  towards  the  surface  but  hitherto 
unreacted  products  are  also  transported  back  in  the  process  in  such  a 
manner  that  they  will  react  closer  to  the  surface. 

X  X  X  X 

The  problem  of  predicting  the  boundary  layer  behaviour  in  a  rocket  motor 
perforation  is  a  very  complex  one  which  does  not  appear  to  he  rigorously 
tractable . 

The  effect  of  mass  addition  alone  is  already  difficult  to  predict:  for 
a  plate,  indeed,  it  results  in  an  increased  thickness  of  the  laminar  sub¬ 
layer  while  for  •  p-»pe  it  may  actually  decrease  the  laminar  sublayer  since 
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the  effect  of  fluid  injection  at  the  wall  is  to  accelerate  the  main 
stream  velocity  f  28  j  .  The  presence  of  severe  temperature  gradients  in 
the  flame  ,  the  complex  grain  design,  and  the  fact  that  the  length  to 
diameter  ratio  of  the  perforation  is  in  general  relatively  small  still 
increase  very  much  the  complexity  of  the  problem. 

Nevertheless,  in  practical  cases  of  erosive  "burning,  the  velocity 
component  perpendicular  to  the  surface  does  not,  generally,  exceed  one 
per  cent  of  the  mean  flow  velocity,  and  it  seems  possible,  following 
Corner  and  Geckler,  to  use  the  well  know  theory  of  flow  in  rough  pipes, 
as  a  first  approximation  [4,  8]  •  Although  this  drastic  simplification 
can  hardly  be  completely  justified  on  theoretical  ground,  it  appears  to 
give  fairly  good  results  in  practice  when  the  theory  is  used  with  caution 
and  experimental  data  are  exploited  for  computing  the  numerical  value  of 
some  parameters  in  an  effort  to  adapt  the  equation  to  the  actual  case 
under  consideration. 

XXX 

In  the  theory  of  flow  in  rough  pipes,  the  mean  flow  velocity  V  is 
expressed  in  terms  of  the  so-called  friction  velocity  v^  by  the  relation¬ 
ship  : 

Rv 

v  =  v _  (1.75  -  5.75  log1Q  — ^  )  (11) 

where  P  is  the  pipe  radius  and  v  =  m /gp  is  the  cinematic  viscosity,  u 
being  the  absolute  viscosity  and  gp  the  specific  weight. 


In  the  laminar  sublayer  the  flow  pattern  is  entirely  determined  by  the 
viscous  stresses.  Neglecting  the  velocity  component  perpendicular  to  the 
surface,  therefore,  it  can  be  written,  by  analogy  with  the  theory  of  pipes* 


y  being  the  so-called  friction  distance  parameter  and  y  the  distance  from 
the  burning  surface.  Equation  (  12  )  defines  the  laminar  sublayer  up  to  a 
critical  value  of  the  friction  parameter  : 


y 


cr 
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which  is  characteristic  of  the  problem  and  represents  the  distance  at 
which  turbulence  appears.  In  the  theoiy  of  rough  pipes,  the  critical 
value  is  about  in  the  problem  under  consideration  it  is  reasonable 
to  expect  a  somewhat  lower  value  due  to  the  mass  addition  and  the  pres¬ 
ence  of  the  flame.  In  fact,  in  a  previous  computation  dealing  with  a 
double  base  propellant  [  91  >  it  has  been  found  t 


and  in  the  following  section,  values  ranging  between  0.4  and  0.7  are 
found,  depending  on  the  ratio  of  perpendicular  velocity  to  mean  flow 
velocity,  for  ammonium  perchlorate  propellants  whose  flame  is  strongly 
heterogeneous  (see  figure  1  8) . 

According  to  the  preceding  assumptions,  therefore,  turbulence  appears 
above  plane  y^  where  temperature  is  reached  in  the  absence  of  sign¬ 
ificant  heat  transfer  from  the  main  diffusion  flame,  and  equation  (12) 
can  be  written: 

*  vt  yi 

XI  t .  -  \ 


Equations  (ll)  and  (13)  provide  the  relation  between  y'  and  V,  assuming 

*  ' 

that  the  proper  value  of  y  is  obtained  from  correlation  of  experimental 

cx 

data.  Elimination  of  v_^  between  these  two  equations  gives  : 


V  =  5.75 


ycrv 


0.30435  +  log. 


£cr' 

1  . 


6  -  Numerical  Application. 


Under  severe  erosion  the  propellant  can  be  schematized  by  a  flat  sur¬ 
face  of  pure  ammonium  perchlorate  whose  solid-phase  decomposition  is 
represented  by  equation  (l  )  or  O’)*  Figure  17  represents  the  surface 
temperature  as  a  function  of  burning  rate,  r,  as  given  by  the  following 
three  recent  references  : 
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A  -  Chaiken  and  Andersen  [  25  1  • 
r  =  4600  exp  (  -  20000/H  Tg) 

B  -  Andersen  and  Chaiken  f  27  1  . 
r  -  31  Tg  exp  (  -  22000/fc  T  ) 

C  -  Johnson  and  Fachbar  t  26  ]  • 
r  -  5.88  Tg  exp  (  -  19160/R  T  ) 

r  being  given  in  cm/sec. 

These  empirical  relationships  have  been  obtained  by  "hot  plate"  pyrolysis. 

Figure  17  also  represents  a  single  point  giving  the  lower  limit  of  T  as 

0 

Friedman  obtained  by  an  entirely  different  technique  consisting  in  the 
optical  measurement  of  the  surface  temperature  [  21  ]  . 

XXX  X 

The  strongly  endothermic  reactions  occuring  at  the  surface  can  be  re¬ 
presented  as  follows  (21,  26]  t 

NH4CI  O^  (s)  -*•  FH^(g)  +  H  Cl  O^Cg)  -  56.OO  kcal/mole  oxidizer  (4) 

For  the  mean  specific  heat  c^  in  the  decomposition  zone,  Johnson  and 
Fachbar  use  [26]  t 

0^  =  38.08  cal/mole  oxidizer,  °K;.  in  agreement  with  our  computation,  at 

1300  °K,  which  gives  =  3t>.34  cal/mole  oxidizer,  °K. 

For  the  present  calculation,  the  following  value  has  been  selected! 

c  =  37.00  cal/mole  oxidizer,  °K. 

P 

Equation  (8),  therefore,  can  be  written  s 

T  -  T'  =  —  =  =  1514  °K 

os—  37 

c 

P 

for  1^=1^=  300  °K,  and  finally  : 

T'  =  -  1214  °K. 
s 

It  is  clear  that  T^  is  not  an  actual  temperature  but  rather  represents 
an  energy  with  its  sign. 

X  X 


X 


X 
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The  adiabatic  combustion  temperature  of  ammonium  perchlorate  is  equal 
to  temperature  at  the  end  of  the  decomposition  zone.  Our  computation, 
assuming  the  simplified  reaction! 

NH4  Cl  04  (s)  -  1.5  HgO  +  H  Cl  +  0.5  N2  +  1.25  0g 
give  s  i 


T1  -  1430  °K 

a  value  which  is  in  agreement  with  » 

T1  >*  1440  °K 

computed  by  Friedman  I  21  ;  at  100  at.  and  taking  into  account  the  minor 
concentrations  of  Cl,  NO  and  NO,,. 

XXX 

Friedman  uses  the  following  value  of  the  thermal  conductivity  in  the 
decomposition  zone  [  21 ] x 

k  =  0.00016  cal/cm  sec  °K 


which  gives  : 


193.7  cm  sec/g 


Our  computation,  using  Eucken's  realtionship,  gives  at  1 300  °K  : 


and 


_£■  •  =  1484  cm  sec/g 

y 


the  following  value  has  been  selected  for  the  present  calculation! 


=  1750  cm  sec  /  g. 


ic 
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Using  these  numerical  values,  equation  ( 1 0) ,  describing  the  heat  trans¬ 
fer  across  the  decomposition  zone,  can  he  written* 

T  +  1214 

-2 -  -  exp  (  -  r  3412.5  y»  )  (lO) 

2644 

for  0^  =  1 .95  g/cm^. 

Used  together  with  equation  (l)  or  (l')>  it  permits  also  to  represent 
thickness  y]  as  a  function  of  r,  for  the  three  solid-phase  decomposition 
laws  which  have  been  considered  previously,  as  shown  on  figure  1,7. 

It  is  seen  that  the  curves  representing  both  T  and  y ' ,  are  relatively 

S  1  -4 

well  grouped  together  *  in  the  following  computation,  curves  B  have  been 
rather  arbitrarily  selected. 


XXX 

Determination  of  the  actual  value  of  the  critical  friction  distance 
parameter  can  be  obtained  through  correlation  of  experimental  data  such 
as  those  of  Lake  and  Marklund  t  14  i  which  can  be  summarized  as  follows* 


Propellant  t  75/25  polysulfide  -  epoxy  t 


p 

cm 

V 

m _ 

sec 

r 

Ofil 

sec 

V 

V 

T 

s 

°K 

7L 

10  H  cm 

100 

150 

1 .257 

0.00703 

1124.0 

0.297 

200 

1 .447 

0.00607 

1140.0 

0.240 

250 

1 .619 

0.00543 

1153.6 

0.200 

50 

150 

0.865 

0.00967 

1082.2 

0.485 

200 

0.965 

000809 

1094.0 

0.420 

250 

1 .053 

0.00706 

1103.5 

0.373 

20 

150 

O.568 

0.01588 

1036.0 

O.85O 

200 

0.603 

0.01264 

1042.5 

0.770 

250 

0.636 

0.01067 

1048.0 

0.720 

where  V 

is  the 

velocity 

component  perpendicular  to  the  surface 

*tli6  rsgi. 

on  of  importance 

for  T  —  1400 

°K  ,  is  given  by  * 

1.5  10 


( cm  /sec  ) 
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981 

8136.6 

27 


h21 

•m3 


r _ 

■ .  P. .  . 

14*00 


8400  r 
P 


P 


being  expressed  in 


a  molecular  weight  of  27  has  been  assumed). 


The  value  of  T  and  y'  have  been  obtained  by  using  figure  17. 

S  1 


The  experiments  of  Lake  and  Marklund  have  been  performed  using  a  tube 
1.5  cm  in  diameter.  In  the  decomposition  zone,  for  T  z  1400  °K,  our 
computation  gives  < 


u 

S  P 


and 


0.562 

981  p  10 


10* 

3 


-3 


0.562 


10 


-3 


2.325 


10-4  p 


g/« 


cm  sec. 


104  1400 


-  2.325  10  4  p 


2.417 

P 


cm 


sec . 


£_ 

3 

cm 


p  being  expressed  in  kg/cm* 


The  value  : 

v  _  2.5  cm^ 

p  sec 


is  used  below,  although  it  is  not  whether  the  mean  value  in  the 

decomposition  zone  or  in  the  mean  flow  has  to  be  used. 


For  T  =  2000  °K  ,  a  rough  estimate  gives  : 

1 0  ^  g/cm  se  c . 

2 

cm 
sec 

It  is  seen  from  equation  (14)  that  y  varies  approximately  as  the 
reciprocal  of  v.  The  selection  of  a  larger  value  would,  therefore,  decrease 
y  accordingly  without  influencing  the  subsequent  computation,  however. 


and 


m  =  0.77 


4.73 

P 
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Equation  (14),  with  the  preceding  numerical  values,  gives  1 

0.30435  +  iog10  i.5yw  +  iog10  jr- 

'  m 

from  which  figure  1.8  has  been  computed,  giving  y  as  a  function  of  V,  /V 

cr  *  * 

for  the  nine  experimental  points  selected  above. 

Figure  18  shows  that  y  seems  to  he  almost  entirely  determined  by 

cr  •  * 

ratio  Vj_  /V,  velocity  V  and  pressure  p  having  only  a  minor  influence.  In 
the  subsequent  computation,  therefore,  the  dotted  mean  line  is  used  for 
approximating  the  dependence  of  y  CT  on  V^/V. 

It  is  interesting  to  remark  that  this  dependence  is  very  much  influenced 
by  the  selected  numerical  value;  in  a  preliminary  calculation  using  : 

Qg  =  36.32  kcal/mole  oxidizer 

a  simple  relationship  between  y  and  V.  /V  was  hardly  noticeable. 

cr 

X  X  X  X 

The  preceding  results,  obtained  by  correlating  experimental  data  with 
theory,  will  be  used  in  the  following  section  for  computing  the  influence 
of  grain  design  and  motor  scale  on  erosive  burning. 


v  y-j  P 

14.375 


cr 
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5.  INFLUENCE  OF  GRAIN  DESIGN  AND  MOTOR  SCALE  . 


Experimental  data  discussed  in  section  3  -  B  suggest  that  local  port 
dimension  covers  "both  the  influence  of  motor  scale  (overall  effect)  and 
that  of  grain  design  (local  effect). 


Quantitatively  the  influence  of  local  perforation  dimension  can  he 
taken  into  account  hy  replacing  the  pipe  radius  R  in  equations  (ll)  and 
(14)  hy  the  local  hydraulic  radius  R^  1 


2  Gross-sectional  area 
Wetted  perimeter 


In  particular,  at  70.3 


kg/c 


cm 


equation  (14)  can  he  written* 


(15) 


V  y» 

0.2045  “  (°« 30435  +  log1Q  ycr  +  log1Q  jt~  )  -  log1Q 

from  which  figure  19  has  been  computed. 


(14) 


It  gives  R^  as  a  function  of  yj  for  twc  particular  values  of  the  mean 
flow  velocity  V  =  200  and  400  m/sec,  taking  into  account  the  dependence 
of  ycr  on  Vx  /V  as  given  hy  figure  18. 

The  dependence  of  burning  rate  r  or  y*  is  also  represented  on  figure 
19  which  serves  as  an  auxiliary  diagram  for  computing  figure  20. 

The  relationship  defining  k^  can  also  he  written  s 


r  -  r 

o 


r 

0 


-  **  °/G" 


(16) 


and  for  a  given  value  of  G/G  (of  V  or  M) ,  in  the  presence  of  a  sufficiently 
strong  erosion,  the  proposed  theory  assumes  that  the  burning  rate  is  in¬ 
dependent  of  rQ  s  it  is  expected,  however,  to  depend  on  hydraulic  radius. 


Comuarison  between  t-wo  cases  differing  hy  strand  burning  rate  r  and 

*  O 

hydraulic  radius  R^  can  he  carried  out  using  the  following  relationship  : 


1000 


FIG. 19 


(17) 
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* 

for  G/G  ,  V  or  V  given 

the  reference  case  being  indicated  by  '  and  the  other  one  by 

In  the  following  computation  no  attempt  is  made  to  arrive  at  absolute 

values  of  the  erosion  coefficient  k,.  but  rather  to  relative  variations, 

M 

assuming  a  reference  point i 

r'  -  0.65  cm/seo 
0 

k'  -  1 .00 

M 

and  R*  «*  1.5  cm. 
h 

at  70.3  kg/cm  .  These  values  are  selected  in  accordance  with  Green's 
data  for  the  purpose  of  obtaining  a  value  of  equal  to  unity: 

The  results  obtained  by  Lake  and  Marklund  suggest  a  larger  value  at 
that  particular  point  while  those  of  Penech  and  Billheimer  may  apparently 
be  interpreted  as  indicating  a  lower  value  (the  motor  dimensions  and  the 
pressure  level  corresponding  to  these  experiments,  however,  are  not  reported) 
t  13,  14,  17J  • 

The  computation  of  figure  20  can  be  carried  out  in  the  following  manner: 

a  -  for  r'Q,  R'^  and  V  given,  the  corresponding  value  of  r'  is  obtained  on 
figure  19  (see  dotted  line'), 

b  -  for  V'  given,  the  values  of  r"  corresponding  to  different  values  of 
R  are  obtained  on  figure  19, 

c  -  for  V  given  and  arbitrarily  selected  values  of  r"  ,  the  corresponding 

values  of  k"  are  computed  through  equation  17. 

M 

Figure  20  shows  that  the  erosion  coefficient  k^  decreases  with  increas¬ 
ing  values  of  hydraulic  radius  IL  and  of  strand  burning  rate  rn. 


r"  (  r'  -  r'  ) 
o  '  o' 

r*  (r"  -  r"  ) 

ft  '  ft  ' 


v* 

r:M 


k" 


M 


p  =  70  kg/cm 


r0*  0.5  cm/ sec 


5  10  20  50 

(cm  ) 


FIG.  20 
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The  range  of  variation  of  k^  is  apparently  in  relatively  good  agreement 
with  the  results  of  Fenech  and  Billheimer,  although  no  definite  conclusions 
can  he  drawn  without  complementary  information  on  pressure  level,  motor 
scale  and  possibly  ballistic  modifiers  [  17  1  • 

For  fixed  values  of  R^,  the  dependence  of  k^  on  rQ  is  normal  i  it 
results  directly  from  the  correlation  for  section  4-C- 6  . 

The  influence  of  mean  flow  velooity  V,  however,  must  be  discussed  with 
the  greatest  caution  )  it  results  mainly  from  the  fact  that  equation  (18) 
is  an  empirioal  approximation  whioh  poorly  describes  the  phenomenon,  and 
does  not  take  into  account  the  existence  in  certain  cases,  of  the  treshold 
velocity  as  discussed  previously  (the  strand  burning  rate  of  the  reference 
point  has  also  been  selected  for  the  fact  that  it  corresponds  to  a 
negligible  threshold  velocity  in  the  interpolation  of  the  results  obtained 
by  Lake  and  Mark}  und  (  1 ,  14  1  )  • 
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6.  A  TENTATIVE  EXPLANATION  OF  IRREGULAR  BURNING  IN  SOLID  PROPELLANT  ROCKETS. 


Several  excellent  survey-papers  have  been  presented  recently  by  McClure  , 
Price  and  others  on  the  important  and  difficult  problem  of  unstable  burning 
in  solid  propellant  rockets  [  29,  30]. 

Generally  these  autho.  s  emphasized  the  poor  present  understanding  of  the 
phenomena  involved  and  the  necessity  to  take  into  account  the  burning 
mechanism  of  the  propellant. 

For  our  purpose,  their  conclusions  can  be  summarized  as  follows: 

-  under  certain  conditions,  many  but  apparently  not  all  propellants  have  the 
intrinsic  ability  to  support  oscillatory  burning, 

-  actual  initiation  of  oscillatory  burning  of  an  intrinsically  unstable 
propellant  results  from  an  extremely  delicate  gain-loss  energy  imbalance 
which  is  influenced  by  a  great  number  of  factors, 

-  in  many  cases  oscillatory  combustion  gives  rise  to  irregular  burning 
characterized  by  large  and  erratic  variations  in  mean  chamber  pressure. 

The  desirability  of  a  clear  distinction  between  oscillatory  and  irregular 
burning  complicates  the  problem  and  the  greatest  caution  is  imperative  in 
this  difficult  matter.  Recent  progi-ess  in  ■understanding  the  steady-state 
burning  mechanism  of  ammonium  perchlorate  propellants  (see  section  4 -A) , 
however,  have  been  sufficient  for  justifying  a  tentative  phenomenological 
explanation  for  the  intrinsic  ability  of  a  propellant  to  support  oscillation 
and  bum  irregularly  under  certain  conditions. 

It  seems  reasonable,  indeed,  to  expect  that  the  surface  erosivity,  and 
possibly  also  acoustic  admittance,  ir.  other  terms  the  propellant  response 
to  small  velocity  and  pressure  perturbations  is  much  larger  for  a  smooth 
surface  than  for  a  rough  one . 

Moreover,  roughness  is  not  the  only  criterion,  another  most  important 
surface  characteristics  being  the  nature  of  the  protrusions.  At  xxw 
pressures,  the  surface  protrusions  are  made  of  ammonium  perchlorate  which, 
alone,  is  a  fair  mor.opropellant  probably  capable  of  amplifying  perturbations. 
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On  the  other  hand,  at  high  pressures,  the  apparent  surface  consists 
in  fuel  peaks  generally  made  of  rather  inert  substances  whose  regression 
can  hardly  he  expected  to  he  influenced  much  hy  small  perturbations. 

This  may  not  he  true  anymore,  however,  for  active  fuels  such  as  nitro¬ 
cellulose  . 

The  discussion  is  summarized  hy  figure  21  which  reproduces  experimental 
data  obtained  by  Landsbaum,  Kuby  and  Spaid  [  31 ]  and  gives  the  severity  of 
irregular  burning,  as  measured  by  the  ratio  of  the  increase  in  pressure  to 
the  normal  chamber  pressure,  as  a  function  of  combustion  to  throat  area 
ratio  K  which  is  itself  a  function  of  pressure  (  K  *»  100,  200  and  300, 
correspond  respectively  to  pQ  =  170,  390  and  710  psi).  This  figure  also 
represents  the  type  of  surface  which  corresponds  to  the  various  behaviours, 
a  nearly  flat  surface  being  assumed  to  give  maximum  instability. 

One  may  be  surprised  by  the  sharp  decrease  in  severity  observed  just 
above  that  particular  pressure.  A  possible  explanation  is  that,  in  this 
region,  the  rapid  appearance  of  even  a  limited  number  of  fuel  peaks  over 
the  main  surface  almost  completely  prevents  the  establishment  of  small 
velocity  oscillations  over  the  active  surface  of  the  oxidizer. 

The  small  amount  of  erosivity  remaining  in  region  D  may  be  due  to  a 
few  large  crystals  still  slightly  protruding  (bimodal  distribution)  and 
in  many  instances  the  existence  of  this  region  has  not  been  observed  [  31  K 

XXX 

This  theory  assumes  that  small  oscillations  mainly  affect  the  protrusions; 
it  is  consistant  with  the  view  that  erosivity  plays  a  predominant  role  in 
irregular  burning  described  by  Brownlee  end  Marble  [32], 

Although  the  acoustic  admittance  of  the  propellant  can  'be  also  expected 
to  depend  on  surface  rugosity,  this  theory  may  not  necessarily  provide  the 
complete  explanation  for  oscillatory  burning,  and  other  factors,  such  as  the 
gas-phase  reactions,  may  have  to  be  taken  into  account,  as  proposed  by 
Barrere  [ 24  1 . 
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Let  us  now  discuss  briefly  some  of  the  conclusions  which  can  he  drawn 
from  this  theory. 

It  can  be  shown  easily,  through  comparison  between  burning  rates  of 
propellants  and  pure  ammonium  perchlorate,  that  an  increase  in  oxidizer 
content  or  a  decrease  in  crystal  average  size  both  result  in  a  decrease 
in  surface  rugosity  and  a  stronger  tendency  for  the  perchlorate  to  stay 
protruding  above  the  surface,  both  factors  favouring  the  ability  to 
sustain  irregular  burning  [ 1 9 »  20]. 

Along  these  lines,  the  influence  of  aluminium  powder  in  suppressing 
irregular  burning  could  be  explained,  tentatively,  by  an  increase  in 
surface  roughness  through  cratering  and  the  detachment  of  crystals  away 
from  the  surface,  as  observed  by  Barrfere  [  24]. 

This  theory  may  perhaps  also  be  applied  to  other  formulations t  indeed 
the  flat  active  combustion  surface  of  double  base  propellants  may  explain 
the  apparent  absence  of  any  particularly  orderly  dependence  of  severity 
of  either  oscillatory  combustion  or  irregular  burning  on  pressure,  in 
tests  on  a  variety  of  formulations  [33]. 

On  the  other  hand  the  systematic  protrusion  of  fuel  peaks  over  the 
burning  surface  of  ammonium  nitre* e  propellants  may  be  the  origin  of 
their  stable  behaviour. 
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